[August 23, 2012] 



Statistics of anomalously localized states at the 
center of band E = in the one-dimensional 
Anderson localization model 

V.E.Kravtsov 1,2 and V.I.Yudson 3 

lr The Abdus Salam International Centre for Theoretical Physics, P.O.B. 586, 34100 
Trieste, Italy. 

2 Landau Institute for Theoretical Physics, 2 Kosygina st., 117940 Moscow, Russia. 
3 Institutc for Spectroscopy, Russian Academy of Sciences, 142190 Troitsk, Moscow, 
Russia 

Abstract. We consider the distribution function P(|?/;| 2 ) of the eigenfunction 
amplitude at the center-of-band (E = 0) anomaly in the one-dimensional tight-binding 
chain with weak uncorrelated on-site disorder (the one-dimensional Anderson model). 
The special emphasis is on the probability of the anomalously localized states (ALS) 
with | if}\ 2 much larger than the inverse typical localization length £q. Using the solution 
to the generating function <I> an (u, </>) found recently in our works |18| I17j we find the 
ALS probability distribution Pd^j 2 ) at |V>| 2 4) 3> 1- As an auxiliary preliminary step 
we found the asymptotic form of the generating function 4> ora (u, 4>) at u 3> 1 which 
can be used to compute other statistical properties at the center-of-band anomaly. We 
show that at moderately large values of I^PA), the probability of ALS at E = is 
smaller than at energies away from the anomaly. However, at very large values of 
\tl>\ £q, the tendency is inverted: it is exponentially easier to create a very strongly 
localized state at E = than at energies away from the anomaly. We also found 
the leading term in the behavior of P(|V>| 2 ) at small \tp\ 2 <C Iq 1 and show that it is 
consistent with the exponential localization corresponding to the Lyapunov exponent 
found earlier by Kappus and Wegner |8] and Derrida and Gardner [3]. 



PACS numbers: 72.15.Rn, 72.70.+m, 72.20.Ht, 73.23.-b 
1. Introduction 

There is a long-lasting interest in localization effects [TJ [2] in Id systems [3]-[18]. The 
simplest and most widely studied model is a linear chain with a nearest-neighbor hopping 
and random site energies e, with no inter-site correlation: (£i£j) = 5ij<j 2 . The wave 
function ipi at a site % of this one-dimensional Anderson localization model [1] obeys the 
equation: 

+ V'i+l + £ i^i = E ^i- (!) 

In the absence of disorder [si = 0) the eigenstates would be plane waves, with 
eigenenergies determined by the wave vector k: E(k) = 2cos(fc), k G (— 7r,7r). In the 
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presence of the disorder, the eigenstates are random and require statistical description. 
Moreover, the states are localized at an arbitrary small disorder strength a. For weak 
disorder the localization length £(E) is large as compared to the lattice constant: 
£{E) ^> 1. This means that the "typical" magnitude of the normalized wave function 
near its localization center can be estimated as \i>\\y P ~ 1/£(E) <C 1. However, for 
some realizations of the disorder, more strongly localized states, " anomalously localized 
states" (ALS), are possible, with the value of the wave function maximum in the range 
of 1/£(E) <C |^| 2 < 1 (the right equality would correspond to a state localized at a 
single lattice site). Our aim in the present paper is to study the probability distribution 
P{\ip\ 2 ) of such strongly localized states in a long weakly disordered chain. 

We will be especially interested in the statistics of ALS in the vicinity of the so- 
called Kappus-Wegner center-of-band (E — 0, k — tt/2 ) anomaly [8J. This anomaly is a 
feature of a discrete chain (it is absent in the continuum model) and originates from the 
commensurability of the de Broglie wavelength and the lattice constant. The anomaly 
manifests itself [HI E] in a sharp, finite in the limit a — > 0, enhancement of the density of 
states (DoS) v[E = 0) and the localization length £{E = 0) inside a very narrow energy 
window (of the width ~ a 2 ) around the band center E = as compared to their values 

u (E = 0) « ±- ; £ = £ (E = 0) = \ (2) 

Z7T (X z 
beyond this interval [19]. In particular, it was shown [9] that in the limit er <C 1: 

u(E = 0) 4V2n 3 £ ext (E = 0) 1 . (\\ 

- 1.01508... ; — y- -/ = -—r 4 - = 1.0942.... (3) 



^o(^^O) T 4 (l/4) ' "' ' £ (E = 0) 16tt 2 U, 
Here we have introduced the superscript u ext" to emphasize that the corresponding 
localization length £ ext = l/[dt^(E)] is defined by the Lyapunov exponent j(E) and 
therefore characterizes the exponentially decaying tails of localized wave functions; for 
this reason it will be referred to as the extrinsic localization length. Similar anomalies 
have been found later [TSJ US] for other physical quantities (like transmission and 
conductance), also related with the Lyapunov exponent. 

In contrast to this set of problems, the eigenfunction statistics P(\ip\ 2 ) may provide 
information about an "intrinsic" spatial structure of localized wave functions including 
the vicinity of the center of localization. In particular, it allows to calculate the 
"intrinsic" localization length £ int (E) = l/h(E), where I 2 (E) = Eil^(^)| 4 is the 
inverse participation ratio. 

However, studying the statistical properties of normalized eigenfunctions is a 
considerably more difficult theoretical problem than studying the Lyapunov exponent 
(the latter is related to propagation of an external wave in a semi-infinite chain and is 
not directly related with eigenfunctions) . 

The formalism for studying the eigenfunction statistics in a disordered chain (see 
review [20]), adapted recently [18J to the case of the center-of-band anomaly, expresses 
moments of the eigenfunction distribution in terms of a "generating function" 0; E) 
of the two auxiliary variables. These variables can be loosely interpreted [H [18] as 
the squared amplitude u ~ |oj| 2 ^o and the "phase" defined by a representation of 
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eigenfunctions in the form: ipj = aj cos (kj + <pj) with slowly varying aj > and 
<f)j G (0, it). The generating function $(u, 0; E) allows one to calculate all local statistics 
of eigenfunctions. In particular, it determines the inverse participation ratio (IPR) 
I 2 and higher moments I m = l]j(|'0j| 2m )> as we U as the full distribution function 
|2 ). Also, the generating function $>(u,(f);E) determines (through a nonlinear 



integral relation Eq.(60)) the joint probability distribution P(u,4>;E) of the amplitude 
and the phase. However, the relationship between the generating function Q(u,<j)]E) 
and the normalized distribution function V(4>; E) = f duP(u, 0; E) of the phase turns 
out to be remarkably simple [18], it is given by the limit u — > of the generating function 
$( M ,0; E): 

<S>{u = 0,<j>;E)=V{<j>;E) = 2P refl {e;E)\ e=2(fl , (4) 

There is also a simple relationship between V(<j)]E) and the probability distribution 
Prefi{9', E) of the reflection phase 9 for a wave incident on a semi- infinite disordered 
chain. It is given by the second equality in Eq.Q. At weak disorder the phase 
distribution V(<f); E) is uniform in the continuum model and outside the center-of-band 
anomaly but it becomes a non-trivial function of at E = [9]. 

A relative simplicity of calculation of such quantities as the Lyaupunov exponent 
(and the extrinsic localization length £ ext (E)) and the DoS, v{E), is due to the fact that 
they can be expressed entirely in terms of the the probability distribution V(4>; E), i.e. 
involve the generating function $(w = 0, 0; E) at u = 0. For instance, the DoS, v(E) is 
given by [HI HH] : 

1; -47t/ #cos 2 (0) [P(frE)] 2 . (5) 



ME) 

On the contrary, the complexity of the problem of local eigenfunction statistics arises 
because it requires the full generating function $(w, 0; E) of the two variables u and 0. 
In particular, the statistics of relatively rare anomalously localized eigenstates of large 
peak amplitude |^| 2 £o ^ 1 which we will study in the present paper is determined by 
$(u, <p] E) at large values of the variable u ^> 1. 

The generating function $(w, 0; E) for a long chain at the center-of-band anomaly 
has been found recently [TTJ dH] by solving exactly the corresponding second order partial 
differential equation Eq.([Io|) in u and <p variables. The exact solution Eq.(14) to this 
equation reflects a hidden symmetry of the problem which has not been yet explicitly 
exploited. However, the solution is given in quadratures as an integral of a product of 
Whittaker functions over the variable which enters both the argument and the index of 
these functions. In this paper we perform a careful analysis of the integral and derive 
the asymptotic form of $(«, 0) at large values of u ^> 1 (from now on we omit the 
energy argument E = for brevity). It has a form: 

Q -^/ub{4>) 

$ an (w,0) = A{4>) ; u>l, (6) 

u i 

where the function 6(0) is a solution to the first order ordinary differential equation 
(54), and A(<p) is specified in the section [2j 
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Figure 1. The logarithm of the ALS probability distribution P(x) (x = \ip\ 2 £o) at 
the center-of-band anomaly E = and outside (E =^ 0). At moderately large values 
of \ip\ 2 lo (H% < 35) the probability of ALS at E = is smaller than that for 
E 7^ 0. However, at IV'PA) > 35 the situation is inverted: the probability (~ 10~ 15 ) of 
very strongly localized states is larger at the band center. Note that at weak disorder 
a <C 1 we consider in this paper the typical localization length £o w 8/cr 2 1 is 
parametrically larger than the lattice constant a = 1 [E5]. Thus the above results 
are valid when the anomalously small localization length is still much larger than the 
lattice constant. 



It allows us to compute the tail of the distribution function P(\ip\ 2 ) at \i/j\ 2 £q 3> 1 
in a long chain of the length L ^> £q: 

^^ ^l^ ' (M 2 4»i), (7) 



where the coefficient k is determined by some "critical angle" C , given by Eqs.(39),(40), 
at which the function 6 2 (0)/4cos 2 reaches its minimum: 

b 2 ((f)r.) 

= 0.830902... < 1. (8) 



4 cos 2 4> c 

The anomalous distribution of eigenfunction amplitudes Eq.([7j) should be compared 
with the "normal" one [18] valid in the continuum model and outside the center-of-band 
anomaly in the discrete chain [23J: 

WWI^ ^y , (H'4»e-*). (9) 

A comparison of Eqs.Q and ^ reveals an unexpected feature (see Fig.l). While the 
probability of moderately strongly localized states (with the peak intensity 1 < |V'| 2 ^o < 
35) is smaller at E — than that away from the anomaly, very strongly localized states 
(with |t/>| 2 A) > 35) are more probable at the band center. Formally this re-entrant 
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behavior is caused by the value of k w 0.83 < 1 (see Eq.(j8j)) at the "critical angle" 
C 7^ 0; for the "normal" case E ^ one obtains 6(0) = 2 and thus C = and k — 1. 

The behavior of moderately strongly localized states is consistent with the result 
Eq. ([3]) for the Lyapunov exponent which gives an enhanced typical extrinsic localization 
length at E = 0. The opposite trend for very strongly localized states is perhaps due 
to the Bragg-mirror effect of the harmonics of the random potential which double the 
period of the lattice [2Tj [22]- 

A point of special interest is the distribution of small amplitudes P(\ip\ 2 ) at 
|VH 2 ^o 1; as it gives an idea on the shape of the tail of the localized wave function. We 
found the leading term |"0|~ 2 in -PQV'I 2 ) at small |"0| 2 ^o an d shown that it is universal 
for all systems with exponentially localized eigenstates. 

The rest of the paper is devoted to the derivation of the announced results and is 
organized in the following way. In section [2] we obtain the asymptotic of the generating 
function $ an (w,0) at « > 1. In the subsequent section [3] we derive the asymptotic of 
the probability distribution function P an {\ip\ 2 ) at \ip\ 2 £o ^> 1- The behavior of -P(IV'I) 
at small IV'I^o is analyzed in Sec. 4. In the last section [5] we summarize and discuss the 
obtained results. 



2. Generating function $ c 



and its asymptotic at u 1 



Sufficiently far from the ends of a long chain, the generating function becomes site 
independent. At the center-of-band anomaly (E = 0) this stationary generating 
function, $ an (w, (ft), obeys the partial differential equation (PDE) [T7| [TB"] 



[1 - cos (40)] u 2 d 2 + sin (40) ud^ + 



3 + cos (40) 







+2 cos (40) ud u — - sin 



2 cos (40) - u 







(10) 



Its solution should also meet the requirements of being a smooth periodic function of 
0, regular, positive and non-zero at u — > (we recall that $(« = 0, 0) is the phase 
distribution function, see Eq.Q) and decaying at u — > 00. 

These requirements are rather restrictive. For instance, the solution 

$o(w, 0) = wexp (— y/u (I cos0| + I sin0|)) (11) 

is not appropriate for it is not a smooth function of 0. 

For comparison, we write down also the equation for the "normal" generating 
function $ norm (w,0) (i.e. when the energy lies outside the anomaly region, or for the 
continuous model): 

3 



u 2 dl -u+~d 2 



± 7 



,(«,0) =0. 



(12) 



This equation looks like a "course-grained" PDE (10) where all the coefficients are 



"averaged" over the angle interval (0, it) (so called "phase randomization") which is 
equivalent to course-graining over the space region £ ^> Ax 3> 1/k. The variables u 
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and in Eq. ( 12 ) are separated and one immediately finds that the only solution decaying 
at u — y oo and remaining regular and non-zero at u — > is given by 



7T 



e~ 2 ^ at u » 1 , 



(13) 



where Ki(x) is the modified Bessel function. This solution has been earlier obtained 
[13] in the continuous model. It also arises in the theory of a multi-channel disordered 
wire [2U [20]. The corresponding phase distribution is uniform: -P norm (0) = <& norm {u = 
O,0) = 1/tt. 

Unlike Eq.(12), the PDE (10) is not separable in the variables u and 0. However, 
due to a hidden (and not well established yet) symmetry of the problem, it was possible 
to find new variables which allowed us to split the PDE ( flo] ) into two ordinary differential 
equation and thus to construct an exact general solution [T7J HH] . The solution, which 
obeys the above requirements, is given by [18J: 

|r(i + eA)| 2 




+ CC. 



(14) 



where e = e l7T ' , e = e~ t7r ' and W\ tfl (z) is the Whittaker function (For the second index 
H = 1/4 the Whittaker function can be expressed also in terms of the parabolic cylinder 
function, see, e.g. |25j). In the limit u — > the expression (14) reproduces the phase 
distribution function P an (0) = $ 6 



u 



0. 



V an {(t>) 



4a/7T 



1 



(15) 



r2 (l) ^3 + cos (40) ' 

which was derived earlier [9] in a different way. It shows that the phase distribution 
becomes non-uniform at the center-of-band anomaly. 

Our current task is to derive an asymptotic expression for $ an (w,0) in the limit of 
large u 3> 1. The integrand in Eq.(14) is too complicated for a brute force attack. This 
is because both the arguments and the first indices of the Whittaker functions are large 
(as is shown below, the leading contribution to the integral comes from A ~ \/u) and the 
standard [22] asymptotic expansions of these functions are not applicable. Our approach 



will include three steps: first we will represent Eq.(14) in the form which allows us to 
find an asymptotic expression of the integrand; then we obtain the asymptotic form 
Eq.Q of the generating function $(«,</>) at large u (this asymptotic expression will be 
obtained in the next subsection), and finally the ALS distribution function P(\ip\ 2 ) will 
be found by a saddle-point integration over 0. 

The generating function Eq.(14) is periodic in (with the period 7r/2) and 
symmetric with respect to the change <fi — > it/2 — 0. Therefore, it is sufficient to calculate 
$(n, 0) in the angular interval (0, vr/4]. We exploit the following integral representation 



Statistics of anomalously localized states at the center of band E = in the one- dimensional Anderson localiz 
of the Whittaker function (cf. 9.222.1 [25]): 



W_, i (x) 



T(l/4 + A) j i 



00 , ff — 1 

e -xt/2 I 1 1 



dt 



t + lj (t 2 -l) 3 / 4 



(16) 



valid for 9ftx > and 9ftA > 0. Since the integrand in Eq.(14) is an analytical function 



within the sector 7r/4 < argA < tt/4, we rotate the integration contour A — > Ae l7r//4 and 
introduce a new integration variable z : 
1 fu 



A 



4 V z 



After these transformations, Eq.(|14|) takes the form: 

2^L 



-3? 



00 e i7r/4 dz 



Here 



where 



fi{t,z,(f>) 



r 4 (i/4) Jo J~z 

dt 



{t 2 - l) 3 /4 



exp [-Vuf m (t,z, 



I ^ /t + 1\ cos 2 . 



4^ \t-l 



is real, while 



f 2 (t,z,4>) = 



1 j + 1 ^ i t \/i sin 2 



4^ Vt-1 



(17) 



;is) 



(19) 



(20) 



(21) 



is purely imaginary for real z. Exact Eqs.(18)-(21 ) constitute the starting point for the 



calculation of asymptotic expressions at u ^> 1. 



2.1. Asymptotic of the integrand in Eq.(18) 



At u ^> 1 the integrals Eq.(19) can be computed in the saddle-point approximation. 



The minimum of the action in the integrand of I\ (z, 0) is achieved at the point 



tr 



1 



z cos z 



> 1. 



(22) 



The integration contour goes through this point, so the corresponding saddle-point 
contribution is given by 



I ( i\z, q 
where 

Fi(z,(j)) = fi(t ,z, 



e -VuFi(z,(f>) 

TiV4 (1 + 2 COS 2 0)V4 

In ( a/1 + z cos 2 + yfz cos 



COS( 

~2 



Z COS( 



(23) 



(24) 



+ y 1 + z cos 2 1 

For the integral ^(z, 0) the situation is more complicated as there are two saddle-points: 

(25) 



t± = ± 



\-r — 1 , z < 1/sin 2 



2 sin (j> 



1 



2 sin iji 



z > 1/ sin 
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Figure 2. Choice of contours (solid lines with arrows) in the complex plane of t: 
a) initial contours; b) and c) the final contours deformed to pass through the saddle 
points. The cuts are denoted by the dotted lines. The contour C_ can be deformed 
away to infinity in the lower half-plane. 



both lie outside the integration semi-axis t > 1. On the complex plane t with two cuts, 
(— oo, — 1) and (1, oo) we define an integral over a contour C by 

r dt 

h[C\ = J c ( t 2_ 1)3 ,4 exp[-V5/ 2 (t,s,0)], (26) 

where we choose the branch of the integrand so that on the upper edge of the cut (1, oo) 
I 2 [C* ] = I 2 (z, 0). Taking the contour C = C L + C\ + C* + as depicted in Fig|2^, 
one checks straightforwardly that 

I 2 [C\ = - 2e^/VvV( 2 v / ^ ([i + ^-vV^y-/^^ 

~ _ 2e^/V^ I/(2 ^ (e l7r/4 I 2 (z, 0)) . (27) 

Thus, with the exponential accuracy we have expressed the quantity of our interest 
3ft (e i7r/4 I 2 (z, 0)) (see Eq.( 

latter is not changed if the integration is extended to parts 6U + and dU^ comprising 
the closed contour (see Figj2^i). In this way we arrive at the important relation: 

h[C] = I 2 [C+] + J 2 [C_] = I 2 [C + ] . (28) 



18 



)) in terms of the contour integral ^[C]- Evidently, the 
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where the contours C+ and C_ are shown in Fig|2b,c; the last equality in Eq.(28) holds 



because the contour C_ may be safely shifted down to the infinitely remote part of the 



half-plane 3t < 0, where -^[C-] vanishes (see Eqs.(26) and (21)). 

Further transformations depend on the location of the saddle-points, i.e. on the 



value of z, see Eqs.(25) and (25). For ,2 sin < 1, the contour C + can be lifted to 



the upper half-plane of t (see Fig[2p) to go through the saddle-point t+ (25) which 



provides the minimum of the action. The corresponding saddle-point contribution to 
3ft (e i7r/4 I 2 (z, 0)) is given by (z sin 2 < 1) 

e -i7r/4 



K(e^ /4 4 s) (z,0) 
where 



-7^/(2^7) As) 



me, 



2tt ^l^g-vAZFj f(z,4 



2mV4 [l 



z sin 



11/4 



sin i 

~2 



1 

z sin< 



7T 1 

— h arctan * / q 

2 V z sin 2 



z sin 



(29) 



(30) 



When zsin > 1, the two saddle points (25) lie on the real axis. We bent the contour 



C + so that it goes through the both points within the proper Stokes sectors (Figj2 
The resulting saddle-point contribution to 3ft (e^^^z, 0)) is given by 

z l/4 



3ft (e iw/4 L 
where 



2tt 
2u W 



smi 



—i 



z smi 



z sin 



In i 



- 1]V4 



3 w/4 -^5^(2,. 



+ c.c. 



z > 



sm 



2; sm 



1 + 



zsmi 



+ 



Z7T 



+ iv z sm 



, (31) 



.(32) 



The two saddle-point expressions for 3ft (e ln ^ 4 I 2 (z, 0)), Eqs.(|29|) and (|3l|), can be 
represented by a single formula valid for an arbitrary z > 0: 



3ft(W 4 4 s) (z,0)) 
where 



27^/4 



2«V4 



3ft 



e i7r/4 e - 



/SF 2 (+) (z 



+ e- 



-17T/4, 



/u-Fi _) (2 



Lz sm 



- 1]V* 



(33) 



F 2 (±) (^, 



sm ( 



In ±i\ z sin 



1 ± sin 



z sini 



z sm 



(34) 



Eqs. (33) and (34) are defined on the complex plane z with a cut along the ray 



(1 / sin 2 0, 00); branches of (z sin 2 — l) 1 / 2 and (z sin 2 — l) 1 / 4 are chosen to be positive 
on the upper edge of the cut, the (standard) branch of In w is defined by the requirement 
S(lnw) = at w > and the cut along (— 00, 0) on the w-plane. Accounting for Eqs. 



(24) and (|33j) , we arrive at the expression for the generating function Eq.(18) in the 



form (u 1): 



7±(Z,0) = 



2ji r Qi-K/4 Q -y/uF+(z,<t>) _|_ e -i7r/ 4 g- 



fuT— (z,i 



\(z cos 2 + l)(^sin 2 0- I)} 1 / 4 



Fi(z, 



(35) 
(36) 
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4.xl(T - 



2.x 10" 



u=800, <j>= 



:/» 



0.002 0.004 0.006 0.008 0.010 0.012 0.014 



Figure 3. The saddle-point integrand (arbitrary units) in Eq.(35) for u — 800 and two 
different angles <p — \ an d 4> — § (thin solid line) compared with the proper integrand 



in the exact solution Eq. ( f 4 ) (thick solid line) . The coincidence is very good except in 
the vicinity of the branch cut point z = sin~ 2 4> where there is an integrable singularity 
in the saddle-point integrand. As u increases, this peak singularity moves to the tails 
of the integrand (to the right tail for <fi < (f> c and to the left tail for <f> > <fi c ) and thus 
makes negligible contribution to the z-integral. An exception is the case of <fi ~ 4> c 
where the peak does not move to the tails. In this case the saddle-point integrand 



Eq. ( 35 ) is no longer valid (see Appendix A) . 



Here the integration contour C = Cq + on the complex plane z with the cuts along 
the rays (— oo, — 1/ cos 2 0) and (1/ sin 2 0, oo), is shown in Figj5] the chosen branch of 
+ 1) 1//4 is positive at z > — 1/ cos 2 0. The integrand in Eq.(35) is depicted (for 



[z cos 



different values of 0) in FigJ3| together with the result of the direct numerical evaluation 



of the integrand (after switching to the z- variable) in Eq.(14). 



Our next step is the calculation of the integral in Eq.(35). 



2.2. Saddle-point calculation of the integral in Eq.(35) for $(«, 



Saddle-points of the integrand in Eq.(35) are determined by solutions z±((j>) to the 
equations 



dJ r ±(z 1 . 

dz 



1 

Iz 



In ( a/ z COS 2 0+1 + y/z COS « 



+ COS (t>\ Z COS 0+1 



±i- 



ln ±i\/z sin 



1 ± i\fz sin 



± % sin0\/z sin 



0. 



(37) 



It turns out that the solutions z±(</>) are real and < z±(<f>) < 1/ sin 2 0; the solution 
Z-((f>) exists for < < C , while the solution z+((f>) exists for C < < vr/4 (we 
recall that we consider the angle interval (0, 7r/4)), see Figj5} The critical angle C is 
determined by the condition z±(0 c ) sin 2 C = 1, i.e. the solution reaches the origin of 
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Figure 4. Solutions z±{(j)) to the saddle-point equations Eq.(37) 



the right cut. With this condition, the equation Eq.(37) results in (for the both signs): 

sin 6 i 



Y((j) c ) = , where Y(<j>) = In 



COS0 7T 



(38) 



v 1 + cos (f) J sin 2 
This transcendental equation can be represented in a nice form using the parametrization 



cot . 



sinh 



(39) 
(40) 



where x c ~ 2.4164... is the solution of the equation 

sinh x — x = it 
The value of the critical angle C 

(t> c = 0.58060... (41) 

arises as an important constant also in the calculation of the probability function 
P an (H 2 ), section H 

In the vicinity of this critical angle we have: 



(42) 



i 



i 



sin (0 C ) sin 



sin 



where Y(<p) is given by Eq.(38). At the ends of the angle interval, i.e. at = 7r/4 and 



0, the solutions to Eq.(37) are given by 

1000 /7T ^ 3 

~ "9" V4 ~ 



z + (cf>) « 20 



7T 



4 



Z-(0) = sinh 2 (xo/2) w 4.1263... ; where sinh(xo) — Xq = 2n 



(43) 
(44) 



Using Eq.(37), one can represent the saddle-point actions in a following compact form: 
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© 



c, 



o 
X- 



-cos 2 o 1- sin 2 



c R 



C R 



a) 



-cos 2 



/ 

Z + 




sin t? 



c R 



C R 



b) 



Figure 5. Contours in the complex plane of z: a) evaluation of the contribution of 
the saddle-point z_ ; b) evaluation of the contribution of the saddle-point z + . At the 
critical angle </> c the saddle points (denoted by a cross) touch the branch cut point 
z = sin -2 4>. 



^±(0) = cos0 \J 1 + z±{4>) cos 2 ± sin0 y 1 - z±{4>) sin 2 . (45) 

In a similar way, the second derivatives of the actions in the saddle point can be 
represented as: 



dz 2 



cos" 



=F 



sin 



=z±(4>) 



cos 2 w 1 — z± (0) sin 



(46) 



It follows from Eq.(46) that the second derivative is positive at 2_(0) (0 < C ) and 
negative at z + ((f>) (0 > <f> c ). Therefore, for < C , the contour C = Co + C? 



sec 



Figj5^i) goes through the saddle point Z-(<f>) within the proper Stokes sectors, and the 
term with T- makes the contribution to the integral in Eq.(35). The term with J r + in 
Eq.(35) does not have a saddle point at < C and its contribution is negligible at large 
u. 

On the contrary, for > C the contour C is not appropriate as it goes within 
improper Stokes sectors of the saddle point z+(<f>). To overcome this obstacle, let us 
modify the contour C = C + C+ by adding an additional contour which corresponds 
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to the lower edge of the cut, see Figj5)D. It is seen easily that this operation does not 
change the integral Eq.(35) because its integrand is purely imaginary on C R . Now, the 
part C R = C+ + C 11 of the modified contour can be deformed to the vertical contour 
which goes through the saddle point z+(<f>) within the proper Stokes sectors, Figj5|D. The 
corresponding saddle-point contribution to the generating function Eq.(35) at > C is 
determined by the J-+(z, 0) term in the integrand. Summarizing these results we arrive 
at the following asymptotic expression for the generating function: 



$:>,0)=A(0)- 



u 



1/4 



U > 1 



Here the function A(<fr) outside of a narrow vicinity of the critical angle 
is given by 



A(0) 



r 4 (i/4) 



sin J 1 + z±(4>) cos 2 T cos 3 J I — z±(0) sin 



1/2 



(47) 
(see below) 

(48) 



while the function 6(0) is given by Eq.(45): 

K0) = n^)- (49) 

In these equations the upper (lower) signs stand for > C (0 < C ). At the particular 
angles, = C , = 7r/4 and = 0, the functions 6(0) and A(<fr) are given by (see 



expressions Eqs.(43) and (|44 
tt/4) 



6(0 
6(0 

A{4> = tt/4) 
A(0 = 0) 

a(<p = C ; 



V2 ; 6(0 = 0) = y/l + MO) 

(50) 2 



2.2641.. 



(50) 



cot . 



823 ^ 



7T2 



r 4 (|; 



+ cot J 



0.6855.. 



2(2vr 



.(0) w 0.3703 



4 a/2 



7T2 



0.6059... . 



sin 



(51) 
(52) 
(53) 



The plots of the functions A((j>) and 6(0) computed from Eqs.(48) and (49) are given in 



Fig|6]and Fig|7j Remarkably, the plots which were calculated from different expressions 
at > C and < C do not show any singularity at = C . The two pieces of the 
curves match perfectly at the critical angle = C . 

In the next section we present a different calculation of the function 6(0) which 
does not possess any critical angle by construction and coincides identically with the 
above saddle-point expressions. As both A(0) and 6(0) are expressed through the same 
solutions z±((f>) of the saddle-point equation, smoothness of 6(0) at = C implies also 
the smoothness of ^4(0). 
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2.3. Ordinary differential equation for the exponent 6(0) 



Let us look for an asymptotic (u ^> 1) solution to the original PDE (10) in the form: 
$(w, 0) ~ exp [— u p b((p)] where m and 6(0) are to be determined by keeping in the 
PDE terms of the leading order in u. We find immediately that p = 1/2 (which is in 
accordance with (47)) while 6(0) obeys the ordinary differential equation (ODE): 

3 + cos40 f db(<j)) \ 2 | sm40 d6(0) | l-cos40 &2 1 

A \ d<j) ) 2 dcj) 4 

One can reduce the equation to the form convenient for numerical integration by 
introducing the function: 

= v^ W 2) 

' 2(l+cos 2 *)i v ' 
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-Log[ c v 



17 

16 : 
15 : 
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13 1 



-Log[<D(u, 0)] 



u=50 



0.2 



0.4 



0.6 



0.8 



given by Eq. ( 47 ) and 



Figure 8. The comparison of the approximate function Q(u, < 
the exact $(it, 4>) Eq.( 14 (obtained numerically at u = 50. Near </> = the exact function 
has an extremum to ensure smoothness of the even function 4>) — $(it, —</>). The 
approximation is not accurate in the vicinity of <fi = 0, where the right branch-cut 
point in Figj5] moves to infinity. However, it becomes more and more accurate as u 
increases. Note that the principle parameter of the approximation it -1 / 4 « 0.38 is not 
very small at u = 50. 



Then Eq.(54) takes the form: 



dy 



-i o sin A 

1 - y 2 



(l+cos 2 <^) 1 /2 



2(1 



cos- 



The initial conditions for Eqs.(54),(56) follow from Eq.(50): 
6(tt/4) = V2, y(n/2) = 1. 



There is an obvious solution to Eq.rt54^ with the initial condition Eq.(57): 



cos d> + sin < 



(56) 



(57) 



(5* 



It corresponds to the choice of sign "+" in Eq.(56). This solution is a growing function 
of (ft with the maximum at = 7r/4. Therefore it does not correspond to the saddle- 
point solution which has a minimum at <\> = it /A (see Fig |7] and Eq.(50)). In fact, the 



solution Eq.(58) corresponds to the particular solution Eq.(Ulh which we have already 



discarded on physical grounds. Thus the relevant solution for our problem is the one 



which corresponds to the sign "minus" in Eq.(56). This ODE can be transformed into 
the Abel's ODE 



but it does not belong to the classes with known solutions. 



We solved Eq.(56) numerically applying the initial condition Eq.(57) at a point 
7r/2 — 5 with 5 = 1CT 10 . We checked that the solution corresponding to the sign 



"plus" matches the function 6q(0) obtained from Eqs.(55),(58) with the same accuracy. 



Much less trivial is that the solution for 6(0) corresponding to the sign "minus" in 



Eq.(56) coincides (with the same accuracy) with the saddle-point solution given by 
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Eqs.(45),(49). Remarkably, the solution to the particular Abel's ODE appeared to be 



represented in terms of the solution to the transcendental saddle-point equation Eq.(37)! 
In this connection we would like to remind about another "miracle" of the problem. 
Namely, the saddle-point solution for 6(0) which we obtained from two pieces J-±((f>) 



expressed through solutions z+(0) and £_(</>) of the saddle-point equations Eq.(37), 
appeared to be smooth at the critical angle C where the two pieces match perfectly. 
Now we understand that this is a direct consequence of the fact that 6(0) can be obtained 
from the ODE which has no singularity at = C . 

This argument is also important to realize that the asymptotic function Q an (u, 0) 



(47) is valid also in the vicinity of the critical angle C where the saddle-point expression 



for the integrand in Eq.(35) is no longer valid. As is shown in Appendix A, at 
|0 — C | < m -1 / 6 the integrand should be modified so that the (fake) singularity at 
z = sin -2 is rounded. Then the z-integral can be computed analytically which results 



in the same asymptotic Eq.(47) with 6(0) given by Eq.(49). We conclude therefore that 
different procedures for |0 — C | ^> u~ x ^ and |0 — C | <C u 1 ^ give the same result. Thus 
there is no real "critical angle" in the function $ an (w,0), while there is a critical point 



in the integrand in Eq.(35). 



3. Probability distribution function Pan{\ip\ 2 ) of anomalously localized 
eigenstates 

The generating function $(u,0;P) allows one to calculate all local statistics of 
eigenfunctions. The probability distribution function P(|-0| 2 ) is connected with a "joint 
probability distribution function" P(w, 0) (see [18J for details): 

p(h 2 )= r°du r^h/f-^os^pWH r ^p(-^fi,A ■ (59) 



o Jo Jo cos 2 • 



The function P(u, 0), in its turn, is related with the generating function $(-u,0). This 
relation in the limit of a long chain of the length L ^> i§ reads: 

P(u, 0) = i^ ^ d u [ +l00+ ° % e 4 ^ $ V, 0) (60) 

Lv{E)U J-ioo+O t 

u (E) r+ioo+o 



Lv(E)u z J-ioo+o 



Lv(E)u 2 J-ioo+o 

where the localization length ^> 1 (away from the E = anomaly), the averaged 
DoS v{E = 0), and the DoS vq(E = 0) of an ideal (without disorder) chain are given 
by Eqs.Q and Our aim is to find the asymptotic form of P{\ip\ 2 ) at IV'P^o ^ !■ 

The asymptotic form of the function P(u, 0) is determined by that of the generating 
function $(m, 0) 

$: n (n,0)=^l(0Ke-v /Sb W ; U >1 (61) 



represented in the form suitable for both the normal (q = 1/4, 6(0) = 2; Eq.(13)) and 
anomalous (q = —1/4; Eq.([47|)) functions. 
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Figure 9. The ratio of the probability distribution P(x) (computed numerically 
from the exact solution Eq.(14|) to the asymptotic expression Eq.(67) obtained for 
x = |V>| 2 ^o S> 1. The depletion at x < 4 arises because of the 1/x behavior of the exact 
distribution at small x as compared to 1/x 3 behavior of the asymptotic expression. 
The quasi-constant behavior at x > 4 (with the value of the constant close to that of 
Eq. ( 68 ) depicted by the dashed line) indicates on the setting up of the pre-exponent 
cx 1/x 3 . 



Plugging this function into Eq.(60) and doing the saddle-point integration over t 
one obtains: 

4g+l/2 



ME = 0) 
Lu(E = 0) 



iq+l u 4<z-3/2 e -ut fe 2 (0)/4_ 



Now using the 7r/2-periodicity of the integrand in Eq.(59) one finally arrives at: 

-4 2 (0) 



P 

where 



C 



C 



2Wg-3/ 2 



tt/2 



4 4 <? L 



ME 



[COS i 

=_0) 



189-1 



4g+l 



We b z (4,) 

4 cos^ (b 



(62) 



(63) 



(64) 



In the limit iV'l^o 3> 1, the major contribution to Eq.(63) comes from the vicinity of 
the minimum of the function 

m\ 2 



B(4>) = 



COS( 



(65) 



entering the exponent. Outside the anomaly (or for the continuum model) the function 
b((p) = 2 (see Eq.(13), so the minimum value of B(<f)) is achieved at = 0. Performing 
the saddle-point integration in Eq.(63) we obtain the announced expression Eq.(|9]) 
for the asymptotic of the "normal" probability distribution function P n orm(\4>\ 2 ) of 
eigenstates. It is interesting that this asymptotic form coincides with the exact function 
for any value of \ip\ 2 - 



Statistics of anomalously localized states at the center of band E = in the one- dimensional Anderson localiz 



For the center-of-band anomaly, the function 6(0) is more complicated. It follows 



from Eq.(50) that the function B(<ft) has its minimum exactly at the critical angle 
(Eqs.(|3§-(|g)) : 

5(0 C + 



1 



+ 



+ ... . 



(66) 

sin" {(p c ) sin' {(p c ) 

Thus the critical angle remarkably appears again in the theory even though we have 
shown that the function Q an (u, (f)) is smooth at (f) = <p c . 



Pn 



where 



The saddle-point integration in Eq.(|63|) leads to the following result: 
1 / 



! ) = C{4> c ) 



C(0 C 



exp 



4 sin 2 



cos 



64tt v / 2 



f4(l) sin 2 c 



3.20. 



and 



1 

4 sin 2 i 



0.8310... < 1 



(67) 



(68) 



(69) 



Equation (|67|) is the main result of our paper. 



4. P(H 2 ) at small \ip\Ho < 1. 

In this section we consider the behavior of the eigenfunction amplitude distribution 
function P{\ip\ 2 ) at small values of l^l 2 ^)- Generically, the small amplitudes \ip\ 2 £$ "C 1 
arise either (i) due to localization when the observation point r in ip = ip{r) lies outside 
the localization volume, or (ii) due to the proximity of the observation point to the 
node of the wave function. In the case (i) the amplitude of exponentially localized 
eigenfunction cannot be smaller than ~ £q 1 / 2 e - L / 2£OXt ; while in the case (ii) the 
amplitude can be arbitrary small. It is clear that the case (i) is realized with almost 
certainty in a large sample, while the case (ii) has small probability proportional to 
the small distance of the observation point from the node. This should lead to the 
drastically different behavior of the distribution function for |^| 2 ^o ^ e~ L//£CXt (case (i)) 
and for \ip\ 2 £o <C e~ L l^ (case (ii)). Our approach based on the exact solution of the 



stationary (with respect to the coordinate along the chain) evolution equation (10) is 
capable of describing only the case (i), as the crossover to the alternative case (ii) and 
the corresponding solution for the generating function are essentially L-dependent. 

Furthermore, one can argue that for the case of pure exponential localization 
the asymptotic behavior of P{\ip\ 2 ) at small \ip\ 2 £o (e~ L / £cxt < \ip\ 2 £o <^ 1) should 
be always P{\ip\ 2 ) = C norm /\tp\ 2 . Indeed, in this case the normalization integral is 
logarithmically divergent and dominated by |^| 2 ^o ~ e - L /^ cx \ Thus the normalization 
constant C norm oc 1/L, as it should be in order to make the first moment (\ip\ 2 ) = \ 
as the eigenfunction normalization requires. Should the profile of the localization tail 
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be of the form L a e L l e> with an extra power-law pre-exponent, the characteristic 
sub-leading terms appear in P(|-0| 2 ): 

Cnorm ,^ Oi 



' norm 
|2~ 



|^| 2 ^ ln|^| 2 + 



1 > 



| 2 £o»e L ^ xt . 



(70) 



Thus studying details of the distribution function -Pd^l) a t small amplitudes one may 
infer information about the profile of the tail of the wave function. 

In this section we briefly discuss how the principle term in Eq.( 70 (arises in our 

2 



at 



< 1 



may arise only when P[u <C 1, 0) oc u . This means that the integral in Eq.(60) 

/■+«oo+0 A+ 

/ fe^V^) (71) 

J— ioo+0 t 

must be proportional to « at « < 1. Then one immediately concludes that the function 
$(n, 0) should have a singularity at u — 0. Indeed, in case of a regular w-expansion, the 
term cx u in $ 2 (u, 
to the integral: 



would result in a linear in u contribution in Eq. ( 71 ) proportional 



+ioo+0 



ioo+0 



dte 4t / e ° = 0. 



To obtain the desired dependence P(u <C 1, 0) oc u 1 , one has to assume that there is 
a term oc ulnu in the expansion of $(u, 0). Then the corresponding integral in Eq.(71) 

r+ioo+0 

(72) 



/-Hoo+u , /-oo , V _ 

dt lnte 4 ^ = -2m / dte-^ e ° = -2m 
-ioo+0 JO 4 



would be non-zero and result in P(u, 0) oc vT 1 . We see, therefore, that the term u lnu in 
the expansion of $(w, 0) at small w is the direct consequence of exponential localization. 
Next, one can check that that the "normal" (away from the anomaly) generation 



function Eq.(13) has, indeed, the expansion with lnw-terms: 

1 oo 



(73) 



(74) 



(U, 0) = - 2 (9n + fn ^ U) U n , 
71 n=0 

where the first few coefficients of expansion are given by: 

go = 1, fo = o, fx = i, si = 27 - 1, 

, 1 5 

/2 = 2' 92=1- 4, 

where 7 = 0.577216... is the Euler constant. 

A natural assumption would be that the generating function at the anomaly 
Q an (u, 0) has the same type of expansion Eq.(73) but with the 0-dependent coefficients 



g n {4>) and / n (0): 



E \9n(<f>) + /n(0) bit] U* 
n=0 



(75) 



If so, one can find the coefficients by plugging the series Eq.(75) directly into Eq.(10). 
Then one obtains the chain of recursive equations: 



3 + cos 40 



(^n — ~ \ sin 40 + n(n — 1) — (n — l)(n — 2) cos- 



#n(0) + 
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+ [sin 409^ + 2n - 1 + (3 - 2n) cos 40] / n (0) = g n -x{(j>) 

d(/, + n(n — 1) — {n — l)(n — 2) cos 40 



3 + cos ' 



-di+ln- 



sin- 



4 

=/n-l(0). 

For n = 1 the equation (77) takes the form 
3 + cos 



fn(<j>) 



-cL — sin 40 



<Vi(0) = o 



(76) 



(77) 



(7? 



which determines the derivative of /i(0): 



(79) 



Witt) = -, 

y/3 + cos (40) 

Now we are going to apply the condition of periodicity of 0) as the function of the 
angle 0. Since the left-hand side of this equation is a derivative of a periodic function, 
its integral over the period must vanish. This can be provided only with the choice 
Cf = 0. Hence, the function /i(0) is a constant: 

fM = F 1 , (80) 



where the constant F\ cannot be fixed by the homogeneous equation (78). 

However, it appears that the requirement of periodicity of the function g\ (0) helps 
to fix the constant F x - Indeed, the equation for this function is: 



3 + cos - 



-d<j> - - sin. 



%i(0) = 9o((f>) - (1 + cos40)F! 



(81) 



where go(4>) = V an (4>) is given by Eq.(15). Looking for the solution in the form 
d(t,gi((f>) = c 9 (0)(?o(0) and taking into account that go((f)) obeys the homogeneous 
equation (for the zero R.H.S.), we obtain the following equation for c g (0): 

4 T 2 (|) 1 + cos (40) 



<Vg(0) 



(82) 



3 + cos 40 ^3 + cos (40) 

Using once again the periodicity condition, we must require the integral over the period 
of each side of the above equation to vanish. This determines uniquely the value of Fx'- 

Fx = ^ . (83) 

We see that the solution can be found uniquely only if one assumes the periodicity (and 
hence smoothness) of the function $(u,0), and this solution corresponds to /i(0) ^ 0. 
This means that it would not be possible to find any periodic solution without a term 
oc u\n.u in the series Eq.(75) for $ an (w,0) . Thus the periodicity requires the singular 
expansion at u = with certainty. We note that the assumption of smoothness was the 
key point to obtain the exact solution Eq.([l4~|) [T7] ITS] . 

With the coefficient fx (0) = F x established one immediately finds the leading term 
in the P{\ip\ 2 ) at small \ip\ 2 : 

T 4 (l/4) i 1 £ ext 1 
16tt 2 L jVf ~ ~L~ |Vf : 



«1). 



4) 
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As was expected, the numerical coefficient in Eq.(84) exactly corresponds to the 
replacement £q — > £ ext in the leading term of expansion of Eq. (|9] in accordance with 
Eq.(§. 

However, finding sub-leading terms in P(\ip\ 2 — > 0) is a separate non-trivial problem. 
We leave its complete study for future publications, outlining here only the origin of the 
difficulties. The point is that one can obtain the formal series of the type Eq.(75) for 
$ an (u,0) with the coefficients g n (4>) and f n (4>) represented by a two-fold integrals. To 
this end we exploit again the integral representation Eq.(16) of the Whittaker functions. 
Plugging it into the exact solution Eq.(14) we do the A-integration exactly using the 
well-known integral 



dX 
A 2 



exp~^~ bA 



£ 

a 



The result is expressed through the two- fold integral: 




dU 



(ti-iy 



55) 



56) 



eti cos 2 



eto sin 2 . 



where e 



,171-/4 p 



eti cos 2 (f) + et 2 sin 



In Eq.(86) one can immediately recognize the combination y/uKi{^/u...) which 
enters Eq.(13) and which generates the series Eq.(75). The coefficients g n {4>) = 

(u. 



are 



fi^O) +g i n ) (<j)) and f n {4>) = g$(<p) (f n /g n ) in the corresponding series for $ fl 
expressed in terms of the coefficients g n and f n appearing in the expansion Eq.(73) of 
2^fuKi(2^fu) and the two- fold integrals: 



9 { n ] (<t>) = CngnRef-^l 



dU 



{4>) = c n f n Re^ — -3^ 



{tl-l 



dU 



{t\-l)*h (tl-l)f 
x L^tiM)[G{ti : h)\ n Wi,h)\ n ~ l 



17) 



(88) 



of the three functions: 
G(ti,t 2 ) 

T^(ti, t 2 



e In 



h + 1 
U - 1 



+ em 



U-l 



where c„ 



eti cos + et 2 sin cp 
L^tiM) = In [G(ti,t 2 ) T^h, t 2 )/8] , 

23(1-") 

r 4 (i/4) • 



(89) 

(90) 
(91) 



One can check that g = g^ ' coincides with the phase distribution function V an ((j)) 
defined in Eq.(15). Furthermore, /i(0) appears to be manifestly ^-independent and 
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coincides with Fi found above (see Eq.(83)). The functions g± (4>) (which is also 

(2) 

independent) and g\ (4>) are also well defined. 



However, starting from n = 2 there is a problem in Eq.(87). As the function 



^(^1,^2) grows linearly with ti i2 , the integrals in Eq.(87) are divergent for all n > 2 



This signals that the expansion Eq.(75) for Q an (u, </>) breaks down. The reason is that 
Eq.(75) does not guarantee the correct, decaying at large u behavior of the generating 
function $(w,0). Thus an additional series in $ an (w,0) may be required to cancel 
possible divergence at u — > 00 of the function obtained by the analytical continuation 



of the series Eq.(75). As the result the sub-leading term in the expansion of $ an (w,0) 
at small u is not proportional to u 2 lnw (as for the generating functions away from 
the E = anomaly) but could be much larger. This anomaly deserves a separate 
investigation. 

5. Discussion and Conclusion 

The goal of this paper was two-fold. The first objective was an asymptotic analysis of 



the exact solution Eq.(14) for the anomalous (at the center-of-band anomaly, E = 0) 



generating function $„„(«, 0) at large values of u ^> 1. The corresponding result is 



expressed by Eqs.(45) and (47)-(49). 



Knowing this asymptotic one can compute various quantities of interest related with 
the local statistics of eigenfunction amplitudes. The simplest one is the distribution 
function of the eigenfunction amplitudes P(\ip\ 2 ) which behavior at large \ip\ 2 gives 
an idea about the probability of anomalously strongly localized states. To find this 
asymptotic form at the center-of-band anomaly was our principal physical objective. We 



managed to obtain the asymptotic expression for P an (\xjj\ ) in a compact form Eqs.(67) 



and ( 68 ) . The result is a bit surprising, as it shows a re-entrant behavior summarized in 
Fig.l, which points out on the two competing physical phenomena behind it. Another 
indication of the same phenomena was first found in our earlier works [T7J [18] where we 
noticed two different scales characterizing the moments of \ip\ 2 - 

We also analyzed the asymptotic of $ an (w,0) at small values of u. The leading 
term $ an (O,0) gives the distribution function of phases V(4>) [HI IS] which is related 
with the distribution of scattering phases. The next-to leading term oc u In u contains 
information about the tail of the typical wave function. We have computed this term 
and shown that it is compatible with the exponential localization with the Lyapunov 
exponent found in Refs. [H [9] . We have also shown how the sub-leading term u In u 
results in the universal leading behavior of P (|"0|) oc \1p\~ 2 at small \ip\ 2 - 

However, it appears that computing the further terms of expansion at small u 
(which contain information about the pre-exponential behavior of the tail of localized 
eigenfunctions) is a non-trivial problem which deserves further investigation. 
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Appendix A. Generating function $ an (w,0) in the vicinity of the critical 
angle <p c 

The saddle-point derivation of the large u asymptotic of the anomalous generating 



function $ ara (-u, 0) in the form (47) with the angular dependent pre-exponential function 



(|6| is valid everywhere except for a close vicinity of the critical angle C , where the 
saddle-point solutions z±(<f)) to Eq.(37) are close to the branching point 1/ sin 2 <f> so that 
the power expansion of the actions J-±(z,(j)) breaks down. 

To find the width of the critical region of small 5<f> = <f> ~ 0c we introduce 
the distances A± between the saddle-points z±((f>) and the right branching point 
z — 1/ sin 2 ((p) : 

A ± = ^ Y --z ± (<P)^(5<P) 2 . (A.l) 
sm m 



The latter estimate follows from Eq.(|42j) in the vicinity of the critical angle; the 
dependence A± ~ (5(f)) 2 is clearly seen in Fig|4j 

The previous saddle-point approach to Eq. p5| ) is justified as long as the width of 
the saddle-point peaks 

\z - z ± (4>)\ ~ u- 1 / 4 \d 2 F ± (z, ( f ) )/dz 2 \- 1/2 ~ ^(Ai) 1 / 4 

is much smaller than the distance Az± from the right branching point. It follows from 
here that \5(p\ should be greater than u^ 1 ^ 6 . 

In the narrow region (of width \8(j)\ < w -1 / 6 ) around the critical angle C , the 
regular series expansion of the actions J r ±(z,<f>) breaks down and the previous saddle- 
point approach is not applicable. To treat this narrow critical region, we will develop 
a modified approach. In fact, we will restrict the analysis to even narrower vicinity of 
4> c : \5(j)\ <C u -1 / 6 , which is sufficient for the calculation of the eigenstates distribution 
function P an (\if)\ 2 ) performed in the section^ 



Note that the integrand in Eq.(35) needs revision in the critical domain, too. This 



is because the saddle point estimate Eq.(33) for the integral I^ZjCp) Eq.(19) does not 



work when the two saddle points t± Eq.(25) approach each other (both go to zero at 
z — y 1/ sin 2 0). The saddle-point estimate of the integral is valid only as long as the 
width of the saddle-point regions \t — t±\ ~ u~ x ^\bz\~ x ^ j sin0 is small as compared to 
the distance \t + — 1_| ~ \5z\ 1 ^ 2 sin between the two saddle points. Here we represented 
the integration variable z in the form: 
1 



sin 2 - 



+ 5z. (A.2) 
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From the above estimates one finds that the saddle-point approach is not applicable 
when \5z\ < u -1 / 3 and therefore \t±\ < vT x l % . For small \b~z\ <C 1 the leading 
contribution to the integral h[C] = h[C+] ( see Eq.(28) and the discussion around 



it) comes from a narrow vicinity of the origin t = 0. Expanding the function f2(t, z, 0) 



Eq.(21) in small t (and keeping only linear terms in 5z) } we arrive at the following 
expression for ^[C]' 



h[C] 



4 l e 



in ip 



dte 



/u- sin q> 



-tSz sin 4> 



(A.3) 



2lZ, 



This expression and the relation Eq.(|27|) determine the quantity of our interest, 



&an(u,(fi) takes the form: 



; the integral representation (18) for the anomalous generating function 



[u,< 



2 1 / 3 (2 7 r) 3 / 2 sin 2 / 3 (0) u 1 ' 12 
r 4 (l/4) ' 



-VS[cot(0)-sin(0) Y(<f>)/2] 



d{8z) e"^ ^ W Y ^ Ai (_ 2 -2/3 sin 8/3 (0) u l/3 Sz 



(A.4) 



with Ai(x) being the Airy function. According to Eq.(38) the function Y(0) ~ 5(p in the 



critical region, so the linear in Sz term in the exponent of the integrand can be safely 
omitted. Indeed, due to the convergence of the integral of the Airy function, an estimate 
for typical 5z in the integral Eq.(A.4) is 5z ~ m -1 / 3 , hence the linear in 5z term in the 
exponent is estimated as u 1 ^ 2 5(p5z ~ v}^8(^ which is negligible in the considered critical 
region \5<p\ w -1 / 6 . Neglecting the linear in 5z term in the exponent and calculating 
the remaining integral of Airy function, we obtain $„„(«, 0) in the vicinity of the critical 
angle C in the form Eq.(47) with the pre-exponential function A((p) — > A c (<f)): 

2(2tt) 3 / 2 



AM) 



(A.5) 



r*(l/4) sin 2 c ' 

This expression matches perfectly the out-of-critical expression Eq.([6]) when the latter 
is formally extended to the critical region (p w C , where z±((f>) ~ l/sin 2 </>. This 
comparison completes our derivation of the asymptotic of the generating function 
$ an (u, 4>) both in and out of the "critical region" and shows that $ an (u, <fi) is a smooth 
function of even in the vicinity of the "critical angle" <p c . 
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